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Abstract. We present a projector quantum Monte Carlo study of the topological properties of 
the valence-bond-solid ground state in the J-Q 3 spin model on the square lattice. The winding 
number is a topological number counting the number of domain walls in the system and is 
a good quantum number in the thermodynamic limit. We study the finite-size behavior and 
obtain the domain wall energy density for a topological nontrivial valence-bond-solid state. 


1. Introduction 

While topological order has mainly been discussed in the context of exotic states such as quantum 
spin liquids mm, similar topological quantum numbers can also arise in systems with long- 
range order. The simplest example of a topological quantum number is the winding number of a 
dimer model 0. In the classical close-packed dimer model on the square lattice (as an example 
of a simple bipartite lattice), a winding number W = (w x ,w y ) can be defined by assigning a 
directionality (arrow) from sublattice A to sublattice B for each dimer. Superimposing any 
such configuration of arrows onto two reference configurations (conventionally one with the 
dimers forming horizontal or vertical columns), closed loops are formed and the winding numbers 
correspond to the total x and y currents normalized by the system length L [3]. 

In the quantum dimer model (QDM), the classical dimer configurations are the basis states. 
For S = 1/2 quantum spins on a bipartite lattice, any total spin singlet can be expanded in 
Valence-bond states (VBs) with positive-defined expansion coefficients. The VBs are dimers with 
the added spin-singlet internal structure and the configurations can again be classified according 
to the winding number |3]. Since the Hamiltonian of the QDM contains only local operators 
that have no matrix element between states of different winding numbers, the Hilbert space can 
be divided into different winding number sectors, i.e., the topological winding number is a good 
quantum number. However, for general VB states, longer bonds are introduced, which makes 
the bond operators in the Hamiltonian non-local. For finite size systems, the winding number 
can be changed if bonds of lengths exceeding 1/4 of the system length appear. Therefore, the 
winding number can only be conserved in the thermodynamic limit. 

To investigate the emergent conserved winding number in an S = 1/2 spin system, we here 
study the J-Q 3 Model on the square lattice |S]- The Hamiltonian of the model is 

H = -J^ Cij - Qs X CijCklCmn (1) 

(ij) (ij,kl,mn) 


with Cij = j — Si ■ S j recognized as a singlet projector. The J term of H is the standard 
Heisenberg exchange and the Q 3 terms have the singlet projectors arranged in columns, as 
shown in Fig. Ufa). This model hosts a columnar valence-bond solid (VBS) ground state when 
Q3/J is sufficiently large 0. Ignoring quantum fluctuations, the VBS ground state is sketched 
as an ordered VBS breaking the Z 4 symmetry of the Hamiltonian, as shown in Fig. Gib). As 
previously studied in the context of QDMs [J], a VBS with domain walls on a periodic lattice 
corresponds to a non-zero winding number—a case of W = (1,0) is shown in Fig. GIc). 
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Figure 1. Illustrations of (a) the Q 3 terms, (b) the columnar VBS state, and (c) a VBS state 
with winding number W = (1,0) (on a periodic lattice). 


In this paper, we consider the pure Q 3 model without any J term to study the VBS 
ground state. We will demonstrate explicitly, using projector quantum Monte Carlo (PQMC) 
calculations, that there is an emergent topological quantum number (a winding number) 
associated with the VBS order in the thermodynamic limit, while in a finite system the winding 
number is not conserved. In this case the topological sectors correspond to different types of 
domain walls, which are unstable (non-conserved) in finite systems but become stable in the 
thermodynamic limit. 

2. Projector Quantum Monte Carlo method 

The model with Q 3 > 0 in (GI does not have sign problems. We thus take advantage of the 
PQMC method with efficient loop update algorithms mmm to study the model. The method 
is based on applying a high power of the Hamiltonian to a trial state IV’o)) 

IV’m) = (-brriV’o). (2) 

For sufficient large rn. I^m) approaches the ground state |0) of the model. As the ground state 
| 0 ) of a bipartite J-Q3 model is guaranteed to be a total spin singlet, it is particularly convenient 
to use a trial state expanded in the VB basis \V r ) in the singlet sector. Then only the singlet 
sector is considered from the outset, leading to faster convergence with m. 

The energy of the projected state can be calculated in the following way m 

(N\ip m ) (N\(-H) m W 0 ) ’ 1 j 

where | N) is one of the Neel states in the 2 -spin basis, which has equal overlap with all VBs. 

By writing {—H) m as a sum over all possible strings of the individual Q3 terms in Eq. 
(HD, the above expression is evaluated by implementing importance sampling of the operator 
strings with the efficient loop algorithm la mg. The implementation of the algorithm to the 
J-Q3 Hamiltonian has been described briefly in |12] . In the end, the energy E is calculated as 
E = —{rid Erif/2) with rif the number of bond flips and rid the number of diagonal operations. 






3. Results 

We chose a basis state \i/jq) = | V r ) as the trial state and consider the cases that \V r ) in different 
winding number sectors W. The detailed definition of the winding number of a VB state can 
be found, e.g., in Refs. mm- In the PQMC simulations, besides the energy, we also sample the 
probability, P(W), of a projected state in the topological sector W = (w x ,w y ). This is done by 
calculating the winding number of each projected VBs Pk\V r ), with P^ a operator string with 
length m generated in the MC processes. 

In the case that the trial state is a VBs in the winding number sector W = (0,0), the 
ground columnar VBS state will be projected out quickly, i.e., within a small m/N , which can 
be defined as a projecting ’’time” (closely related to imaginary time P3]). This is indicated by 
the convergence of the ground state energy density eo (L) for a system with linear size L. 

Now turn to the cases that the trial state is in the nontrivial winding number sector 
W 7 ^ (0,0). For small systems, the columnar VBS state is again projected out after some 
projection time m/N. This is indicated by the convergence of the energy density e(L) to the 
ground state value eo (L). Meanwhile, the probability P(W) decreases (to 0 for L —> oo). 
However, as the system size increases, the projection time m/N needs to grow as well in order 
for the ground state to be obtained. In the thermodynamic limit, we expect that the system 
will stay in the sector with the initial winding number W, and then it is also plausible that the 
energy of the system will converges to a value ew > eo corresponding to the lowest excited state 
within the sector W. 

To demonstrate such behavior, we introduce the energy density ew(L) of states in the winding 
number sector W, which is obtained by only sampling those states in the sector W. Figure [2] 
shows the ’’time evolution” of the probability P(W ) for the system staying in the original winding 
number sector W = (0,0), W = (1,0), and W = (2,0) as a function of m/T 2 , respectively, for a 
system with linear size L = 96 (lower panel). The corresponding energy density ew(L) converges 
to the values which are higher than eo (L), if W / 0 (upper panel). 



Figure 2. The probability P(W) and 
the energy density ew{L) as functions of 
time m/N (projector power rescaled by 
the system volume). 



Figure 3. Snapshots of {B a ( r)) for a periodic 
system with winding number W = (1,0), in which 
a 27 t domain wall (four separate tt/2 domain walls) 
is formed (upper panel) and for an open system 
with appropriate boundary conditions in which a n 
domain wall is forced (lower panel). 


We now study the reason of the energy gap between a system in a nontrivial topological 
sector and in the W = 0 ground state. It is well known that the ground state of the J-Qz is the 
columnar VBS. The VBS state can be detected by the columnar VBS order parameter, which 

















































































is defined by the operators 


= B £ (x, y)(-ir, Dy = ±Yl B v^ y^ V - 

x,y x,y 

where the dimer operator B a is used: 


B a { r) = S(r) • S(r + a) 


(4) 

(5) 


where a = x,y denotes the lattice unit vector in the two directions. In a columnar VBS, either 
D x or D y has a nonzero expectation value. 

However, if a state bears a nontrivial winding number, the lattice translational symmetry 
with period two is also broken. As an example, Fig. [3] (upper panel) shows a snapshot of the 
bond configuration (B a (r)) for a projected state with winding number W = (1,0). To describe 
such a state, it is useful to define the local order parameters for the VBS state. Consider the 
case that there is still translational invariant VBS with period two along one axis of the square 
lattice, say y axis. The one-dimensional local x and y order parameters as a function of the x 
coordinate are defined m as follows: 

D x (x) = [{B £ (x,y)) - ^(B £ (x - l,y)) - ^(B £ (x + l,y))](-l) x , (6) 

D y (x) = [{B$(x,y)) - (B${x,y + l))](-lf. (7) 


These two local order parameters are independent of the y coordinate. Following theoretical 
expectations uni, The VBS angle 6 (x) can also be defined |I3] , 

Dy(x ) + Dy(x + 1 ) 

2 D x (x) 

such that 0 = 0 and 0 = n for a fully x or y oriented VBS order, respectively. 

To understand the energy gap between system in an nontrivial topological sector and in the 
W = 0 ground state, we first set asymmetric open boundary conditions on a L x L cylindrical 
lattice, where the x-direction boundaries have been modified. By removing Q 3 terms with 
vertical bonds closest to the left edge on every even row and closest to the right edge on every 
odd row, we obtain a state with forced VBS angle 0(0) = n/2 and 0(L) = — n/2. A snapshot of 
the bond strength (B a ( r)) for the state projected out from simple trial state with a corresponding 
domain wall is shown in Fig. 0 (lower panel). More properties of the projected state are shown in 
Fig. EH (right), where the two local order parameters and the VBS angle are plotted as functions 
of x coordinate and various system sizes L. In this case, the local order parameter D x (x ) is 
maximized in the center and goes to zero at boundaries; D y (x) is about 0.75 at the left boundary, 
but gradually changes to —0.75 at the right boundary. The VBS angle 0(x) changes from n/2 
to —n/2. Clearly a domain wall is formed, as expected. According to the total change of the 
VBS angle, we define this as a n domain wall (which consist of two n/2 domain walls). 

The energy density of the state with n domain wall (L) can be calculated as described 
above. The resulting e n (L) is larger than the ground state energy density eo (L) of a L x L 
cylindrical lattice with symmetric boundary modification (i.e., forcing a VBS with no domain 
wall). The energy difference is entirely due to the presence of the domain wall. We thus define 
the domain wall energy per unit length as 


0(x) = atan 


K n (L) = [e n (L) - e 0 (L)\L (9) 

As shown in Fig. O when system size goes to infinity, the domain wall energy per unit length 
converge to a constant K n , which is estimated as 0.434(9), as listed in Tabled) 
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Figure 4. The local order parameters D x (x ), D y (x ) and the VBS angle 6 {x) for projected 
states in the topological sector W = (1,0) with periodic boundaries (left) and for the states 
with a 7 r domain wall (right) in a system with asymmetric open boundaries. 


Similarly, a 7r/2 domain wall can realized by setting asymmetric open boundary conditions 
forcing horizontal and vertical bond patterns at the edges. The energy difference to eo (L) 
multiplied by the system size again converges to a constant when system size L tends to infinity 
(see Fig. ED. The constant is thus identified as the domain wall energy n ^/2 f° r the 7 t/2 domain 
wall with the estimated value being a half of the 7r domain wall energy, as listed in Tab. |T| 

We now return to systems with periodic boundaries. In the case of large enough system, the 
projected energy eigenstate is still in the winding number sector W of the initial state. Take the 
W = (1, 0) case as an example. Figure |4] (left panel) shows the local order parameters and VBS 
angel as functions of the x coordinate. D x (x ) has a minimum in the center and maximums at 
the left and right sides, while D y (x ) is almost zero at the edges gradually grows to a maximum 
of 0.5, then decays to a minimum about —0.5, and finally returns to 0. As the result, the VBS 
angle 0(x) changes from 0 to 7r/2, then to — 7r/2, and back to 0. Domain walls are clearly 
formed here and their dynamics in a large system is very slow and practically locked in to their 
locations in the initial state. According to the total change of the VBS angle we define this as 
a 27r domain wall (which can be regarded as four elementary 7r/2 domain walls). A finite-size 
scaling analysis of the energy density gap to the ground state (W = 0) ew{L) — eo (L) multiplied 
by the system size L is shown in Fig. El Clearly, the scaled gap converges to a constant, which 
can be understood as the 2-k domain wall energy per unit length K 2 n- The estimated K 2 n is 
indeed twice that of the 7r domain wall, as seen in Table [T] 

We further calculated the energy density gap ew{L) — eo(L) from the projected state in the 
winding number sector W = (2, 0) to the ground state for several system sizes. The local order 
parameters and VBS angles as functions of the x-direction shows that there is a Air domain wall. 
Again the energy gap multiplied by the system size converges to a constant, which is the 47r 
domain wall energy. The estimated domain wall energy is 0.205(9), which is about four times 
of the 7r domain wall energy. 
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Figure 5. Finite-size scaling analysis of the domain wall energies. The curves are fitted 
polynomials, used to extrapolate the energies listed in Table [Q 


domain wall 

tt/2 

7 r 

2vr 

47r 

m/N 

6 

6 

4 

3 

K 

0.0208(5) 

0.0434(9) 

0.093(9) 

0.205(9) 


Table 1. Domain wall energy per unit length for various domain walls. m/N is the projection 
“time” at which the energy is calculated in each case. 


4. Summary and Conclusions 

Using a ground-state projector QMC method in the VB basis, we have studied projected states 
in various winding number sectors of the J-Q 3 model on the square lattice. We showed that 
the projected state stays in the winding number sector of the trial state when system size goes 
to infinity. The energy of the states in a nontrivial winding number sector have gaps to the 
VBS ground state in the W = (0,0) sector, due to the presence of domain walls (the winding 
number counting these domain walls). Such domain walls are stable only for infinite size, thus, 
the winding number is an emergent quantum number in the thermodynamic limit. 
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